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ABSTRACT 


Thin liquid films when perturbed, either go to rupture or obtain an equilibrium 
configuration due to interfacial instabilities combined with evaporation (or condensa- 
tion). This behaviour is studied using a continuum approach ( Navier-Stokes eqation 
modified for intermolecular interactions). The dynamics is studied by numerical sim- 
ulations of the non-linear partial differential evolution equation. Born repulsion is 
also incorporated to study the dynamics of the film even after rupture. The process 
of formation of nano-drops starting with a uniform thin film is studied. Usually hy- 
drodynamic instabilities evolve on a short time scale, which first lead to a (quasi) 
mechanically stable structure. Evaporation becomes significant at longer times. 

This technique is further applied to study the generation of critical nucleii under 
heterogeneous nucleation and the energy associated with these. An attempt is made to 
create the critical nucleii having minimum energy associted with them. Since critical 
nucleii are very small (~ 1 nm), the finite size effects including excess intermolecular 
interactions cannot be ignored as is done in classical theory of heterogeneous nucle- 
ation. Our results show that the energy barrier for nucleation may be many orders 
smaller than the predictions of the classical theory, which implies much faster rates 
of nucleation. 
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Chapter 1 

Introduction And Literature 
Survey 


Thin liquid films are those which are thin enough so that their thermodynamic and 
transport behaviour depend on the film thickness. These films are in the range of 
10 - 1000 A. The interfacial hydrodynamics of these films subjected to various kinds 
of instabilities (long and short range intermolecular interactions, Marangoni flow, 
thermo capillarity effects. Rayleigh- Taylor instabilities etc.) is governed by all these 
forces. Their stability analysis gives the kinetics of deformation and equilibrium 
configurations. 

Thin films are encountered in a lot of industrial applications, e.g. rupture of 
foams and soap bubbles, coalescence of drops or emulsions, enhanced oil recovery, 
froth floatation, wetting of surfaces, heterogeneous nucleation and in biomedical ap- 
plications (e.g. cell adhesion to substrates, tear film rupture etc.). 

Nucleation comes into picture in the formation of a new phase [1]. Small clusters 
of molecules form and these grow by accretion to the point of becoming recognizable 
droplets or crystals that may finally coalesce or grow to give massive amounts of new 
phase. For very pure substances nucleation does not take place easily because of the 
extra surface energy of small clusters. In the presence of a foreign material (e.g. a 
solid surface), the tendency of a liquid having contact angle less than 180° is to wet 
the surface. This wettability of the liquid on the solid makes it easier to attatch to 
the solid and ,therfore. a smaller nucleus can grow. It reduces the energy required to 
create the critical nucleus. 

The linear theory for falling films was first studied by Yih [2] and Benjamin [3]. 
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Benny [4] and Atherton and Homsy [5] derived the non-linear evolution eqation at 
the film interface. Scheludko [6] and Vrij [7] studied the effect of van der Waals inter- 
actions and estimated a critical film thickness at which a non-thinning film becomes 
unstable and ruptures. Ruckenstein and Jain [8] modelled the liquid as a Navier 
Stokes continuum modified by van der Waals interactions. They used the lubrication 
approximation to obtain the dynamic linear stability results. Williams and Davis 
[9] adapted the long wave reduction procedure and derived the non-linear partial 
differential equation describing the evolution of film interface. Teletzke, Davis and 
Scriven [10] put forward a theory of wetting of solids by liquids accounting for capil- 
lary pressure gradient, gravitational potential gradient, surface tension gradient and 
disjoining pressure gradient. Stability of free films was studied by Sharma, Kishore 
and Salaniwal [11]. 

The effect of interfacial mass transport on the stability and flow of thin filmswas 
studied by Wayner [12] and an expression for interfacial mass flux was obtained. 
Similar expressions were obtained by Sharma [13]. For evaporating films vapour recoil 
instabilitis were studied by Hickman [14] and Palmer and Maheshri [15]. Bankoff 
[16] studied the evaporating film assuming the interface to be always at saturation 
temperature. Burelbach, Bankoff and Davis [17] considered an evaporating film with 
heat transfer and derived a model for evaporation. They predicted complete dryout 
(film rupture) for positive Hamaker constant. Sharma and Ruckenstein [18] employed 
a non-linear analytical theory which was in agreement with the results of Burelbach 
et al. [17]. Analytical findings of Sharma and Ruckenstein [19]- [20] show discrepancies 
in rupture times as computed from linear and non-linear theories. The effect of acid- 
base forces in polar systems has also been studied [21]. The role of intermolecular 
interactions in coexistence of drops with thin films was studied by Sharma [22] and 
Sharma and Jameel [23]. 

In the first part of the present work, Bom repulsion is incorporated to the inter- 
molecular interactions to study the films even after rupture. In the second part this 
approach is applied to study the generation of critical nucleii under heterogeneous 
nucleation and energy associated with these. 
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Chapter 2 

Evaporating-Condensing Thin 
Liquid Films 


2.1 Introduction 

The evaporation from or condensation to a thin liquid film depends not only on the 
relative saturation in the vapour atmosphere but also on hydrodynamic conditions of 
the film viz. thickness of the film and its intermolecular interactions with the surface 
supporting it. In this chapter we have shown the effect of these instabilities on the 
evolution of the film. The effect of Born repulsion is also studied. The instability due 
to vapour recoil was found to exert negligible effect on the behaviour of the film [18]. 
Therefore, these are not considered. 


2.2 Formulation Of Problem 

2.2.1 Dimensional Equations 

The time evolution eqn. (non-linear hydrodynamic eqn. neglecting the vapour recoil 
instabilities) for a thin liquid film (evaporating or condensing) resting on a plane 
surface (fig. 1) is obtained from Navier-Stokes equation [24] neglecting gravitational 
effects and considering intermolecular interactions. The equations are solved in two 
dimensions only assuming the film to be infinite in the third dimension. The equation 
is given [18], [13] as- 


, m 1 
ht + - - r- 
P %P 


h 3 ( P + <f>) x 


( 2 . 1 ) 


w 


.-here p is the pressure as given by Young- Laplace equation for capillarity- 
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( 2 . 2 ) 


p-p„ = 


“ cf h x 


(i + hif 12 

4> is the force due to intermolecular interactions given as- 


= -2 S LW ^ 
‘ dh h 3 


5 P (do - h' 
exp 


(2.3) 


A — V A ) 

S LW and S p are the apolar and polar components of the spreading coefficeint, S 
[25], [26]. The mass loss term, m, is given as- 


m = a [p 0 exp (j£ (- <?h xx + </>)) - p„] (2.4) 

Here in eq.2.1 , the second term corresponds to the mass loss, the first term in 
the square braces corresponds to the capillary forces and the second term gives the 
intermolecular interactions. The expression for intermolecular interactions takes into 
account only Lifshitz-van der Waals forces and acid-base forces. We have also included 
Born repulsion which we shall see later in this chapter. 


2.2.2 Non-dimens ionalisation 


The evolution equation above contains many physical constants of film liquid such as 
p, p, S LW , S p etc. These can be combined together. The non-dimensionalisation 
is done to study the effect of large number of independent dimensional parame- 
ters in terms of less number of independent non-dimensional parameters. The non- 


dimensionalisation of the evolution equation is done using the following scales- 


H-- 


X- 


Hq 

(1 AjCa) 1 / 2 

h 0 


d ~h 0 


1= 


A 


ho 

^A^Ca^jt 


The following nondimensional numbers are also introduced- 


3piA 


-6 S LW d 2 


Ca= ^K °Ca 

m and cj) are non-dimensionalised as- 


5^ 

JLW 


1 

6 d?P 


1 m f 2 tRT\ 2 

p 0 a \ M u . ) 


i $h 0 
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The evolution equation in its non-dimensional form is obtained as- 


Ht + C\M + H 3 ( Hxxx ~ sgn (A) = 0 


(2.5) 


M — exp [C 2 ( —sgn(A)Hxx + $)] — Ri (2.6) 


where, 


$ = \h~ 3 -PI exp 

O 


d-H 

l 


The non-dimensional parameters obtained are- 


C\ — a 



C 2 = 


(ZApv 2 



(2.7) 


(2.S) 


(2.9) 


Ri = - (2.10) 

P» 

The term sgn(A) implies the sign of A. The surface tension term tends to stabi- 
lize the film whereas the L-W and acid-base components of the spreading coefficient 
[i.e.S LW and S p ) describing the intermolecular interactions exert stabilizing or desta- 
bilizing effect depending on their sign. 


2.2.3 Inclusion of Born Repulsion 

There are three types of intermolecular interactions that come into picture for thin 
films. These are the L-W forces, acid-base forces and the Born repulsion. L-W forces 
and acid-base forces have already been taken into account. These show a rupture of 
the film if the energy change at the minimum distance of approach, d 0 is negative. 
But in reality at this distance the Born repulsion always comes into picture which 
makes the energy change positive below this distance. Here we discuss the criterion 
for inclusion of Born repulsion, and the changes to be made in the expression for the 
intermolecular interactions to include the Born repulsion. 


• When to include Born repulsion : 
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The expression for <p without Born repulsion is given by eq.2.3 as- 


A O cLW^o S P f d 0 — h\ 

* = Sh = - 2S y ~ t p \~ir~ ) 

If the energy change is becoming more and more negative as we come to d Q from 
infinity, it will lead to rupture of the thin liquid film ( if Born repulsion is not ac- 
counted for). Here 4> is positive at this point. And, therefore, it becomes the critereon 
for deciding about the inclusion of Born repulsion. 

So we have to include Born repulsion if 6 ( d 0 ) is positive, i.e., 


.o^lw ( d 0 -d 0 

dl A P I A 


> 0 


i.e., 


S LW ^ ck H q 

< 9A l lf S 


P : 


is 


negative) 


S p ' 2A 

We also have to include Born repulsion if any of the S LW and S p is zero and the 
other is negative (in this situation the A G is always negative). 


How to include Born Repulsion : 


The Born repulsion is included in intermolecular interactions to satisfy the follow- 
ing two constraints at the bare surface (the surface is assumed to be bare when the 
thickness of the film is 1.37 A)- 


A G(l 0 ) = S LW + S p 

(2.11) 

and 0(/ o ) = 0 

(2.12) 

This is done by solving the above equation separately for the polar and apolar 
forces. The energy for intermolecular interactions is written as- 

A G = A G lw + A G p 


For polar forces including Born repulsion, 


A G P =c»*p( A ) + 4 s 

(2.13) 
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where Ci and ci are constants to be determined with the following constraints- 

A G(l 0 ) = ( 2 . 14 ) 

and <j>(l 0 ) = 0 

For apolar forces including Born repulsion. 

AG- = | + |i ( 2 . 15 ) 

where C3 and c 4 are constants to be determined with the following constraints- 

A G{1 0 ) = S LW ( 2 . 16 ) 


and (p{l 0 ) = 0 

The solution of these equations gives the coefficients as- 


Cl 


1 - 5 


-1 


1 

8A 


c 2 = -S p r~ I 1 - A- 


lo 

8A 


-1 


03 = jS™ll 


ci = -5 S™Il 

These coefficients are further modified by multiplying these by another set of 
coefficients <21,02,03 and a 4 respectively so that we don’t incorporate Born repulsion 
for positive values of S LW or S p . These coefficients are each equal to 1 if both S LW 
and S p are negative. If S p is positive and S LW is negative, then 

^ = 1 -—, o 2 =0 

If S LW is positive and S p is negative, then 
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a 4 = 0 


“ 3 = 4’ 

The equation for <j> after incorporating Born repulsion and replacing l a in terms 
of d 0 ( 4 l 3 * * = 3 <P 0 ) is obtained as- 

* = (l - _1 exp + 


_81_ aiS p V0J5g£ ( V0l5d o ' 

256 A h 9 \ 8 A 

Now we nondimensionalise <f> by- 


-l 


2 a 3 S 


LW<P ° +^-a 4 S LW<Po 


h 3 32 


h 9 


(2.17) 


$ 


4>h 0 


6S LW d 1 2 

and the nondimensional form of cj> is written as- 


$(H) = PI 1 


V(L75d' 

3 / , 


-1 r 


(V0l5d-H\ 81VT75 d 9 

-ai exp f 7 


1 a 3 9 a 4 d 6 * * 

3 W ~ 64 If 9 " 


(2.18) 


where P is non-dimensional number obtained while non-dimensionalising the evo- 
lution equation. This is the final expression for nondimensional $ which we use in 
eq. 2.5 and eq. 2.6 instead of eq. 2.7. The evolution eqaution is discretised in a set 
of ordinary differential equations (O.D.E.’s) and then solved. This is discussed next. 


• Discretisation of Evolution Equation 

The initial condition used to solve the evolution equation is- 


H(X, 0) = l + €cos(fcZ) 

The periodicity interval of the governing eqation 2.5 which is from 0 to A is converted 

to the periodicity interval 0 to 2i r by defining a new axial co-ordinate 77 , such that 

77 = kX , where k is the wavenumber of the perturbation. 


8 



The profile of the interface can be defined by any number of nodes and the evo- 
lution eqn. at ’j’th node can be rewritten as- 


Ht{]) - -CiM(j) - k H 3 (j) (k 3 H vvr ,(j) - k.sgn(A).$ r ,(j)') 
The second term on R.H.S. can be rewritten as- 


k 2 


H 3 {j) (k 2 Hr,r,(j) - sgn(A).$(jf) 


Fll(j) - F12(j) 


At] 


where 


Fll(;) = H 3 U + i) (l J n m (j + 5 ) - sgn( A).m + I) 


F 120 ') = H 3 U - 5) - 5) - *gnU).*lj - 5 ) 

Now these are further discretised as- 


F11(j) = H 3 (j- 1 -) 


F21(J_± 1)-F21(j) 

A;? 


where F2hj) = k 2 H m (j) - sgn(A)${j) 


and 


F12(j) = H 3 (j - i) 


F 22(j) - F 22(j - 1) 

A T) 


where F22<j) = k 2 H nri (j) - sgn(A)$(j) 

The second derivative H m is calculated as- 

H(j + l)-2H(j) + H(j- 1) 
vrt ~ (A r ,) 2 

The evolution eqution is now fully discretised where everything at any point of 
time is a function of the height of the film at various points. 
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2.3 Linear Stability Analysis 

The base state is taken as an evaporating (or condensing) flat film [21]. Denoting the 
base state quantities by an overbar, the equation for base state (from eq. 2.5) is- 


Hf + Ci exp (C 2 $) — 


(2.19) 


Linearizing eq.2.5 around base state for vanishingly small amplitude of disturbance 
H' such that- 


H = H(T) + H' 

the equation for time evolution of the perturbation H' will be given by [21]- 

^ + {¥h + sgn(A)k 2 } [c^ + sgn{A)k 2 H 3 ] = 0 (2.20) 

where, T = exp (C 2 $) 

For a disturbance of the form- 


H' = exp (u>T + ikX) 

where uj = growth coefficient, k = non-dimensional wave number, 
the growth rate (coefficient) will be given by- 


u) + + sgn(A)k 2 } + sgn(A)k 2 H 3 


= 0 


( 2 . 21 ) 


where H = H 0 . 

The dominant wave number, k m> at which surface deformations grow at maximum 
rate ( du/dk = 0) is given as- 


k 2 = 


CiC 2 % + H 3 $ h 


( 2 - 22 ) 

2 sgn(A)H 

Comparing eq.2.22 with dominant wave number derived by Sharma and Jameel[25], 
it is found that the film considered here for water at 20°C prefer negligibly larger 
wavelength deformations largely due to mass loss. 
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2.4 Results and Discussion 


• Numerical Schemes Used : 

The non-dimensional evolution equation which is dicretised in O.D.E.s is solved 
using the D02EJF routine of the NAG Libraries. This routine solves the O.D.E.s by 
Gear’s method. 

2.4.1 Equilibrium Film Thickness 

It has been shown earlier [21] that a thin film undergoing evaporation attains a time 
invariant flat profile due to the balance of attractive (or repulsive) intermolecular 
forces and mass gain tor loss). The condition for equilibrium film thickness is (m = 0) 
which together with eq.2.4 gives- 


, ter p v 
© = jf - In — 
Vl Po 


This equation can be rewritten in the form- 


S LW \ (dl\ 1 


S p 


where 


p 


d, Q h G 


A 


= c 


1 RT p v 

c s p V L n 


It is clear from the above equation that there are multiple roots of this equation 
as it is non-linear in h.. Figs. (2) and (5) show these roots for different degrees of 
saturation. From fig. 2) it is seen that for c ~ 0.01 or less, there are even three roots. 
If S LW is positive and S p is negative, then the highest and the lowest of these are 
stable equilibrium thicknesses and middle one is unstable. This is shown in fig. (3) 
that a flat film attain; a stable equilibrium by evaporation or condensation. The case 
studied is for the value of c = 0.0001 and -S LW / S p = 2.0. S LW = 15.0 dynes/cm and 
S p = -7.5 dynes/cm.. If the initial thickness of the film is above the unstable root, it 
will come to the stable root on the upper side through evaporation or condensation. 
Also, if the initial thickness is less than the unstable root, it will attain the stable 
equilibrium on the lower side. If the signs of the values of S LW and S p are reversed, 
the behaviour of the roots is also reversed. This is shown in fig. (4) where the both 
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the cases are compared. For negative S LW and positive S p , the highest and the lowest 
root are unstable while the middle root is stable. Again from fig. (5) it can be seen 
that for negative values of ’c’ from 0 to -1.0, there are two different roots. The case 
studied is —S LW /S P = 2.0, S LW = 15.0 dynes/cm and S p = -7.5 dynes/cm, and c 
= -0.001. For these values, the lower root is stable while the higher root is unstable, 
fig- (6)- If the sign of the values of S LW and S F is changed for the same value of 
c, the upper root becomes stable while the lower root becomes unstable, fig. (6). If 
both the values of S LW and S p are negative, their is only one root which is unstable, 
fig. (6). 

2.4.2 Instability of Water Films ( S LW > 0, S p < 0) 

The non-dimensional time evolution equation is solved using the finite difference tech- 
nique. In the previous work [21] Fourier collocation technique was used to discretise 
the nonlinear partial differential equation. It was shown that a film of water of mean 
film thickness of 4 nm at dominating wave number when perturbed slightly undergoes 
rupture. In the F.C. technique there are some non-physical oscillations at the largely 
flat portions which cause a spurious rupture. Therefore the finite difference technique 
has been employed to study these cases. This technique is found to be successful as it 
clearly shows the ” morphological phase separation” of the film (’’morphological phase 
separation” means the formation of a drop in equilibrium with a thin film). This is 
shown in figs. (7) and (8). The same case is also studied with mass loss. It only 
hastens the growth of perturbation. The drop is formed alongwith a thin film and 
due to mass loss the drop decays to a flat film of eq. thickness. This is shown in figs. 
(9) and (10). 

2.4.3 Effect of inclusion of Born Repulsion 

Those cases in which the spreading coefficient is negative, the tendency of a perturbed 
film is to grow in the form of a drop, leaving bare surface on other areas. But at lower 
thickness, Born repulsion comes into picture (as already discussed when to include 
Born repulsion). It does not allow the thickness to go below the thickness at the bare 
surface and thus the evolution of drop and the hole formed can be further studied 
even after a bare surface is observed. 
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The cosine wave perturbation of Methylene Iodide film is studied. Here S p = 0 
and S LW = —ve. Here we have to include the Born repulsion. Figs (11) and (12) show 
the formation of a drop which further decays to give the bare surface. The maximum 
angle of the drop with horizontal increases very fast as the drop forms and then it 
decreases with the decay of the drop. This is shown in Table (1). 

The same case is studied at higher degrees of saturation (p v /p 0 = 0.8). Due to 
higher saturation, mass loss is less. Hence, the drop formed is of larger size and it 
also takes more time to decay to give the bare surface. This is shown in figs. (13) and 
(14). Again it is seen that the maximum angle of the drop with horizontal increases 
very fast as the drop forms and then it decreases with the decay of the drop. This is 
shown in Table (2). 
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Fig. 2.1 Physical configuration of a thin film lying on a flat surface 
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Fig. 2.7 Film profiles at different intervals of non-dimensional time for a thick 
film of water neglecting mass loss with h 0 = 4 nm. (S LW = +15.46 
dynes/cm, S p = -30.58 dynes/cm, p v /p 0 — 0.5, Ci — 53.576, C 2 = 
-2.6807 x 10~ 4 , P = 9.3907 x 10 3 , k = k m = 2.696(30.0519^m- 1 )) 
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Fig. 2.8 Non-dimensional time evolution for maximum and minimum film thick- 
ness for a thick film of water of h 0 = 4 nm neglecting mass loss. (S LW 
= +15.46 dynes/cm, S p = -30.58 dynes/cm, p v /p 0 = 0.5, C x = 53.576, 
C 2 = -2.6807 x 10- 4 , P = 9.3907 x 10 3 , k = k m = 2.696(30.051 9pm- 1 )) 
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Fig. 2.9 Film profiles at different interv; 

film of water with. h 0 = 4 nm. 
30.58 dynes/cm, p v /p 0 = 0.5, 
P = 9.3907 x 10 3 , k = k m = 2.e 





Fig. 2.10 Non-dimensional time evolution for maximum and minimum film thick 
ness for a thick film of water of h 0 = 4 nm. (S LW = +15.46 dynes/cm 
S p = -30.58 dynes/cm, p v / Po = 0.5, C x = 53.576, C 2 = -2.6807 x 10~ 4 
P = 9.39 07 x 10 3 , k = k m = 2.696(30.0519/im -1 )) 
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Fig. 2.13 Film profiles at different intervals of non-dimensional time for a thin 
film of methylene iodide with h 0 = 2 nm. (S LW = -20.0 dynes/cm, S p 
= 0 dynes/cm, p v /p 0 = 0.8, Ci = 0.2945, C 2 = 1.2384 x 10~ 2 , P = 0, 
k = k m = 0.7058(42.8466^m- 1 )) 
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ig. 2.14 Non-dimensional time evolution for maximum and minimum film thick 
ness for a thin film of methylene iodide of h 0 = 2 nm. (S LW = 
20.0 dynes/cm, S p = 0 dynes/cm, p v /p 0 = 0.8, Ci = 0.2945, C 2 = 
1.2384 x 10~ 2 , P = 0, k = = 0.7058(42.8466/im~ 1 )) 
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Table 2.1: Evolution of a film of Methylene Iodide (h 0 = 2 nm, S LW = -20 dynes/cm, 
S p = 0, p v /p 0 = 0.5) 


Time 
(non- dim.) 

ax 

(non-dim.) 

H m in 

(non-dim.) 

@max 

(degrees) 

Volume 

(nm 3 /nm) 

0.00 

1.1000 

0.90000 

0.6956 

206.9927 

0.40 

1.0491 

0.82330 

0.7844 

194.6849 

1.68 

0.9133 

0.32460 

4.1546 

155.0271 

1.76 

1.7811 

0.06784 

41.1801 

152.5958 

1.84 

3.6574 

0.06784 

42.6166 

151.1661 

1.92 

3.7069 

0.06784 

43.6565 

150.2107 

4.00 

3.3167 

0.06784 

42.4440 

126.2408 

8.00 

2.6574 

0.06784 

42.3941 

84.5866 

12.00 

1.8406 

0.06784 

41.2241 

47.5116 

16.00 

0.7501 

0.06784 

37.6601 

18.9100 

1 20.00 

0.06784 

0.06784 

0.0000 

14.0407 




Table 2.2: Evolution of a film of Methylene Iodide (h 0 = 2 nm, S LW = -20 dynes/cm, 
S p = 0, Pv/Po = 0.8) 


Time 

(non-dim.) 

H ma x 

(non- dim.) 

Hm in 

(non-dim.) 

@max 

(degrees) 

Volume 

(nm 3 /nm) 

0.00 

1.1000 

0.90000 

0.6956 

206.9927 

1.60 

1.0557 

0.70548 

1.2704 

186.9852 

2.40 

1.0516 

0.40356 

3.5844 

176.7264 

2.60 

3.9741 

0.06822 

43.1171 

174.4047 

4.40 

3.8901 

0.06822 

43.7611 

164.1590 

4.80 

3.8360 

0.06822 

43.3039 

161.9120 

13.00 

3.2071 

0.06822 

42.9510 

117.0160 

33.00 

1.2503 

0.06822 

41.8087 

28.0880 

37.00 

0.06822 

0.06822 

0.0000 

14.1200 




Chapter 3 

Heterogeneous Nucleation : 
Circular Critical Nucleii 

3.1 Introduction 

This chapter essentially deals with the situation where the critical nucleus is assumed 
to be a circular drop lying on a flat film. The study is made by considering the initial 
profile of a drop lying on a flat film and then see whether it grows or decays. By 
this study we are able to find a critical nucleus which will grow under the conditions 
prevailing . 

3.2 Generation of Initial Profile 

The initial profile is taken as a part of a circular drop lying on a flat film of thickness 
1.37 A (fig. 1). The equilibrium contact angle 6 is found from the relation : 

cos (6) = 1 H — 

7 

where 5 is the spreading coefficient (S = S LW + [25], [26]. 

Different initial profiles are generated for different heights of the drop (h 0 ) above 


the film. 



3.3 Formulation Of Problem 

The time evulution eqn. describing the behaviour of the drop is the same as shown 
in previous chapter. The only difference is that there is no need to consider the wave 
number as we are dealing with a different initial profile now. The non-dimensionalisation 
is done w.r.t. the top thickness of the drop above the liquid film, h 0 . The evolution 
eqn. is : 


H t + CiM + [# 3 (. Hxxx - sgn (A) § h H x )\ x = 0 (3.1) 

and 

M = exp [C 2 (-sgn (A) H xx + $)] - R\ (3.2) 

The initial profile is defined as in the fig. (1). The growth and decay of different 
nucleii with varying size are studied. The results are discussed in the next section . 

3.4 Results and Discussion 

To find a critical nucleus we study the growth/decay of a drop. The top thickness 
of a drop (h 0 ) is varied. This gives different nucleii of varying radius of curvature. 
A typical case study is shown for the parameters : S LW = 15 dynes/cm, S p = -30 
dynes/cm, p v /p 0 = 1.0224. The study shows that the nucleii with radius of curvature 
greater than 315 Awill grow and those with radius of curvature less than 292 Awill 
decay. This is shown in fig. (2) and (3). It is observed that the drop initially grows 
both in vertical and lateral dimensions. When it grows sufficiently in the lateral 
dimension, it starts interacting with the neighbouring drops and grows as a flat film. 
If we consider the drop to be lying in a larger pool of flat film as shown in fig.(4), the 
growth of the drop in vertical as well as lateral dimensions is more before it starts 
interacting with the neighbouring drops and grow as a flat film. This is evident from 
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comparing fig.(2) and fig. (4). 

The study is also made for two more cases with S p = - 40 dynes/cm and -70 
dynes/ cm while other parameters remaining the same. In the former case the critical 
nucleus is found to have a radius of curvature between 276 AAand 285 A, while in 
the latter case the critical nucleus is found to have a radius of curvature between 262 
Aand 284 A. The equilibrium contact angle increases with S p becoming more -ve , 
which means a larger critical nucleus. This is as predicted by the simulation results 
and shown in fig. (5). 
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Chapter 4 


Heterogeneous Nucleation : 
Energy Considerations 

4.1 Introduction 

In the previous chapter we have seen how to obtain a critical nucleus which is of the 
shape of a circular drop. But so far no consideration has been given to the energy 
associated with it. The critical nucleus need not be necessarily of a circular shape. 
It may have any shape. Critical nucleii of different shapes have different amounts of 
energy associated with them. It is, therefore, of great importance to find a critical 
nucleus which requires minimum amount of energy for its generation. This chapter 
deals with the generation of a critical nucleus with minimum energy, observing the 
growth or decay of drops near the critical nucleus and study the effects of various 
parameters on the size of the critical nucleus. 

4.2 Formulation of the Problem 

4.2.1 Generation of the Critical Profile with Minimum En- 
ergy 

When the mass loss or mass gain at every point of a profile is zero, it means that it 
is an equilibrium profile. This may be a stable or unstable equilibrium. An unstable 
equilibrium is one which when slightly perturbed, either grows or decays depending 



on the perturbation. The mass-loss term at any point is given by: 


. _ ( m w y/ 2 r 

m ~ a \2^RTj [ p ° exp 

At equilibrium, m = 0, which gives: 


Vl ( —crh x 


RT \(l + /i 2 )3/2 


+ <t> 


Pv 


(4.1) 


h xx (1 + h 2 x ) 3/2 + t = o 

v 'LI 7 


(4.2) 


As a flat film of thickness h e is also a solution of this equation, and h. xx = 0 at h e , 

so we get- 


(j>(h e ) RT , p v 

— t — = y 7 n = ( a cons ^ an ^) (4-3) 

Tlu; solution of this equation gives the eq. thickness, h e . 

The equation4.2 when integrated, gives us 

(l + hl)~ 1/2 + ^-K 0 h + c = 0 (4.4) 

or h x = {{K 0 h - - c)- 2 - l] 1 / 2 (4.5) 

7 

This equation is used to generate the critical profile, details of which will be 

discussed while explaining the results. Now, to obtain c, h 0 and R 0 , we proceed as 

follows: 

At x = x c (where the drop is at its maximum height, fig. 1), we have h = 
h 0 and h x = 0. Therefore, from equation4.4 

l + ^M-K 0 h 0 + c = 0 (4.6) 

7 

At x = -oo, we have h = h e and h x = 0. Therefore, from equation4.4 

1 + Ag (M _ KoK + c = o (4.7) 

7 
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From these two equations, we get: 


AG(h e ) + jK 0 (h 0 — h e ) = A G(h a ) (4.8) 

lhe solution of this equation gives the top thickness of the drop, i.e., h a Also we 
get the constant c from any of the two equations above. 

At x = 0, we have h = h 0 , h x = 0 and h xx = =± 

Jrt 0 

therefore, from equation4.2, we get: 


-1 

R 0 


4>(ho) 

7 


+ K 0 = 0 


or 


R 0 = ( K 0 — 


MM \-i 

7 


So we get the radius of curvature of the critical nucleus at the top. 


(4.9) 


4.2.2 Calculation of Energy Required to Generate the Crit- 
ical Nucleus 


The energy required to create a nucleus is given by the difference of the energy 
associated with a nucleus lying on a thin film (fig. 1) and that with a flat film only. 
For a flat film: 


Fi = (71 + 7») A + AG (h.) A - AhM In ^ (4.10) 

y L Vo 

For the drop lying on a flat film: 


F e = 2 7 , [ C {/l + hl)dx + A'Jsl + 2 A G{h)dx - 2^- In ^ ° hdx (4.11) 

The energy required to generate the nucleus is then. 


A F = F C - Fi 


AC\ 



or 


A F 


2 J 0 AG ( h )dx- AG(h e )A + 2 7 ; + hl)dx - 2 1l x c 


2 RT 
V L 




hdx — h e x c 


or 


A F 


/•* C ] r , 

2 y o (AG(h) - AG (h e ))dx + 2li J q (yl + h^-l )dx 


2RT p v 
~rf In 
y l Po 


J rx c 

1 ( h — h e )dx 

o 


( 4 . 12 ) 


where the integral / 0 * e (yl + H%dx gives the area of half the drop. 

This is the energy required to generate a profile over a flat film, the profile being 
defined from x = 0 to x = x c . This equation is used to give the energy required to 
generate critical nucleii, both from continuumm approach and the classical approach. 


Continuum Approach 

The critical nucleus which requires minimum energy for its generation is given by 
the solution of equation4.5 in the range x = 0 to x c using ‘d02ejf ’ subroutine from 
tlieNAG librararies (to be discussed in the results section). Equation4.12 is then used 

to get the energy required as: 


A F cr « 


2 I*' AG{h)dx - AG (h e ) A 
. Jo 


+ 


H jf C {yjl + hl)dx - 2j l x c 


2 RT. p v 
— = — in — 
V L Po 


‘ rx c 

/ hdx — h e x c 
.Jo 


(4.13) 
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Classical Approach 


In the classical theory, the drop is assumed to be large enough so that we can consider 
the AG (h) to be negligible (which is not true for smaller drops as we will see in the 
results). The drop is assumed to be a circular one (fig. 2). a: c is replaced by R x .r c 
is the radius of curvature of the circular drop, # is the equilibrium contact angle and 
AG(h e ) is replaced by 7; (cos 9 — 1). So from equation4.12 we get the energy required 
to generate a critical profile of circular shape as: 


A F = -\r 2 [2 9 - sin 29} £ln^ + 2 r sin 9 
1 V L Po 


-yi9 


sin# 


— 7; cos 9 


For A F to be minimum, = 0, which gives radius of curvature as: 


r , = 


2 H si nfl(sb — c°sa) 

(20 -sin 20) fling 

and the corresponding energy required is: 


(4.14) 


AF* = -ir* [2# - sin 2#] 2 - In V - + 2r c sin 9 

4 Vl Po 


sin# 


— 7/ cos # 


(4.15) 


The energy required for a critical nucleus as from continuum approach and classical 
approach are then compared which is discussed next in section Results and Discussion. 


4.3 Results and Discussion 


• Numerical Schemes Used to Solve Different Equations. 

Different subroutines from the NAG Libraries have been used to solve different 
equations. In addition to the D02EJF routine for solving stiff O.D.E.s, two more 
routines are used - D01GAF and C05ADF. 

. D01GAF 


AO 



This routine integrates a function which is specified numerically at 4 or more 
points, over the whole range specified. It uses third order finite difference formulae 
according to a method due to Gill and Miller. 

• C05ADF 

This routine locates a zero of a function in a given interval by a combination of 
the methods of linear interpolation, extrapolation and bi-section. 

• Effect of Equilibrium Contact Angle : 

The increase in the equilibrium contact angle, #, means that the wettability of 
the surface is decreasing. This reduces the effect of surface on the heterogeneous 
nucleation and, therefore, we need a larger critical nucleus which can grow. As the 
size of the critical nucleus increases, it supports our assumption in classical theory 
that A G (h) can be neglected. So the energy required as predicted by the classical 
theory) AA*) comes closer to the that predicted by the continuum approach. A F cr n 
(shown by simulation results). 

We have studied the effect of changing different parameters on the energy required 
to generate the critical nucleus both from the classical theory and from the continuum 
approach. The study is made for three different situations: 

1. Liquid - Water , S p = 0 , S LW = -ve ; 

2. Liquid - Octane , S p = 0 , S LW = -ve ; 

3. Liquid - Water . S p = -ve , S LW = +ve ; 

The study of all these situations supports the concept that any change that leads 
to an increase in 9 results in larger critical nucleus and also less discrepancy between 
the energy required to generate the critical nucleus from both the theories. 


1. Liquid - Water , S p = 0 , S LW = -ve 



When S p = 0 and S LW is negative, then the spreading coefficient is negative 
and we have to include Born Repulsion. As it is evident from the relation cos (0) = 
l 4 1 1 } and AG (h e ) = S LW + S p in this case, that as S LW becomes more and 

more negative, the contact angle 9 increases which results in larger and larger critical 
nucleus. This concept is clear from fig. ( 3 ) in which we have shown the critical 
nucleus for different values of S LW . 

Next we study the effect of increasing p v /p 0 (higher saturation). At higher satura- 
tion a smallei critical nucleus is required which could grow. This is again supported 
by the results as shown in fig. ( 4 ) and ( 5 ). Different cases have been studied for p v /p 0 
ranging from 1.1 to 2.0. 

The discrepancy in the energy required to generate a critical nucleus as predicted 
from classical theory and from continuumm approach decreases at higher 9 and vice- 
versa. The discrepancy is measured in terms of the ratio of the A F CTit and A F*. 
We have seen that as p v /p 0 increases, the size of the critical nucleus decreases which 
does not support the assumption of classical theory that A G(h) can be considered 
negligible. Therefore, the discrepancy increases with increase in relative saturation. 
Also, when S LW is made more and more negative, the contact angle increases which 
implies a larger critical nucleus size and, hence, less discrepancy. These concepts are 
supported by the results shown in fig. (6) . 

Next we show how a slight perturbation in the critical nucleus affects its growth 
or decay. The behaviour of the nucleus is studied in the same way as in previous 
chapter except that the non-dimensionalisation is done w.r.t. the top thickness of the 
drop above the flat solid surface. Fig. ( 7 ) shows the growth of a nucleus the height 
of which is 5% more than the critical nucleus at every point. Initially the drop grows 
and then it becomes a flat film (due to interactions with neighbouring drops). Then 
it grows as a flat film. Fig. (8) shows the decay of a nucleus the height of which is 
1% less than the critical nucleus at every point. The drop continuously decays and 



becomes a flat film of equilibrium thickness, h e . 

The results of someof the parameter variations are shown in table 1. 

• 2. Liquid - Octane , S p = 0 , S LW = -ve 

Another case is studied where the liquid is Octane. Here also we have to include 
Born repulsion. The variation in contact angle is studied and shown directly in fig 
(9). It is showing clearly that the discrepancy decreases as contact angle increases. 
Also at the same contact angle, if surface tension is increased , S^ w becomes more 
-ve correspondingly ( to keep 0 same) . This results in the reduction in wettability of 
surface. Therefore the critical nucleus size increases and discrepancy decreases. Two 
different sets of curves are drawn for 'p v /p 0 = 2 and p v /p 0 = 3 . The discrepancy 
increases with increase in p v /p 0 for the same reasons as cited before . Fig. (10) shows 
the discrepancy over a wider range of p v jp 0 . The increase in the size of critical nucleus 
with S LW becoming more and more negative is shown in fig.(ll). Figs. (12) and (13) 
show the growth and decay of a critical nucleus if it is slightly perturbed. The growth 
is observed for 5% higher perturbation at every point and decay is observed for 1% 
lower perturbation at every point. 

The results of some of the parameter variations are shown in table 2. 

• 3. Liquid - Water , S p = -ve , S LW = +ve 

When the effect of different parameters is studied in case of water taking polar 
forces into account, we don’t have to include Born repulsion in the parameters range 
studied. Fig.(14) and (15) show the decrease in discrepancy with increase in the 
equilibrium contact angle. Also, the discrepancy increases with increase in relative 
saturation. Fig.(16) shows the same thing for a higher range of p v /p 0 . Fig (17) and 
(18) show the growth and decay of a critical nucleus if it is slightly perturbed. The 
growth is observed for 5% higher perturbation at every point and decay is observed 



for 1% lower perturbation at every point. Results of some of the parameter variations 
are shown in table 3. 



Fig. 4.1 Physical configuration of nucleus of liquid on flat solid surface. 
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X = Xc 
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Fig. 4.2 Physical configuration of a circular drop lying on thin film on a solid 
surface. 






Fig. 4.3 Shape 
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Fig. 4.5 Shape of critical micleii for water for different values of p v /p 0 . (S p = 
0, S LW = -15 dynes/cm) 
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Fig. 4.10 Ratio of energy for critical nucleii of Octane from simulation and clas- 
sical theory for higher range of p v /p 0 . (S p = 0, S LW = -8 dynes/cm.) 



















Fig. 4.14 Ratio of energy for critical nucleii of water from simulation and classical 
theory vs. equilibrium contact angle. (S LW = +15 dynes/cm) 










Fig. 4.17 Growth of a nucleus of water above critical nucleus. (S p = -50 dynes /cm,' 
S LW = +15 dynes/cm, p v /p 0 = 1.2, Cx = 1.5970, C 2 = -0.3794x 10“ 2 , 

P = 444.07) 
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Table 4.1: Results for critical nucleii of Water (S LW = -15 dynes/cm, S p = 0 ). 


He 

Po 

r c 

(A) 

^Gclass 

(ergs/cm) 

K 

(A) 

K 

(A) 

^■G critical 

(ergs/ cm) 

^max 

(degrees) 

Eggs*! 

1.1 

56.59 

0.703 x 10" 4 

1.379 

12.858 

0.606 x 10" 5 

29.26 

0.86 

1.2 

29.58 

0.367 x 10" 5 

1.390 

7.179 

0.268 x 10~ 5 

25.10 

0.73 

1.3 

20.56 

0.254 x 10" 5 

1.400 

5.233 

0.156 x 10 -5 

21.84 

0.61 

1.4 

16.03 

0.197 x 10" 5 

1.412 

4.231 

0.101 x 10~ 5 

19.08 

0.51 

1.5 

13.30 

0.163 x 10~ 5 

1.423 

3.609 

0.681 x 10- 6 

16.64 

0.42 

L8 

11.48 

0.139 x 10" 5 

1.436 

3.180 

0.470 x 10~ 6 

14.42 

0.34 

1.7 

10.16 

0.122 x 10" 5 

1.449 

2.861 

0.325 x 10“ 6 

12.36 

0.27 

1.8 

9.18 

0.109 x 10~ 5 

1.464 

2.610 

0.221 x 10" 6 

10.40 

0.20 

1.9 

8.40 

0.987 x 10" 6 

1.480 

2.403 

0.145 x 10~ 6 

8.52 

0.15 

2.0 

7.78 

0.898 x 10~ 6 

1.499 

2.225 

0.890 x lO" 7 

6.66 

0.10 

2.25 

6.65 

0.705 x 10~ 6 

1.577 

1.816 

0.708 x lO" 8 

1.66 

0.01 

2.28 

6.54 

0.675 x 10" 6 

1.599 

1.748 

0.228 x 10" 8 

0.86 

0.0034 
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Table 4.2: Results for critical nucleii of Octane. ( S LW = -8 dynes/cm, S p = 0 ) 


Pv 

Po 

r c 

(A) 

AG class 

(ergs/cm) 

h e 

(A) 

Hq 

(A) 

AG critical 

(ergs/cm) 

$ max 

(degrees) 

&G*clas3 

1.1 

152.55 

0.131 x 10~ 4 

1.370 

57.309 

0.128 x 10~ 4 

90.00 

0.97 

1.2 

79.75 

0.687 x 10~ 5 

1.372 

30.557 

0.653 x 10~ 5 

44.36 

0.95 

1.3 

55.42 

0.477 x 10- 5 

1.374 

21.599 

0.443 x 10~ 5 

42.73 

0.93 

1.4 

43.21 

0.372 x 10~ 5 

1.376 

17.094 

0.337 x 10- 5 

41.38 

0.90 

1.5 

35.86 

0.309 x 10 _s 

1.378 

14.373 

0.273 x 10~ 5 

40.22 

0.88 

1.7 

27.40 

0.236 x 10~ 5 

1.381 

11.231 

0.199 x 10~ 5 

38.30 

0.85 

2.0 

20.98 

0.180 x 10" 5 

1.385 

8.827 

0.144 x 10“ 5 

36.03 

0.80 

2.5 

15.87 

0.136 x 10~ 5 

1.391 

6.897 

0.991 x 10~ 6 

33.21 

0.73 

3.0 

13.23 

0.114 x 10~ 5 

1.396 

5.888 

0.762 x 10“ 6 

31.08 

0.67 

3.5 

11.61 

0.994 x lO' 6 

1.400 

5.258 

0.621 x 10~ 6 

29.37 

0.62 

4.0 

10.49 

0.897 x lO" 6 

1.404 

4.821 

0.525 x 10~ 6 

27.93 

0.58 

6.0 

8.11 

0.690 x 10- 6 

1.418 

3.874 

0.323 x 10~ 6 

23.79 

0.47 

9.0 

6.62 

0.558 x 10~ 6 

1.433 

3.254 

0.200 x 10" 6 

19.85 

0.36 

12.0 

5.85 

0.490 x 10- 6 

1.446 

2.924 

0.140 x 10" 6 

17.11 

0.29 

15.0 

5.37 

0.446 x 10~ 6 

1.457 

2.709 

0.104 x 10- 6 

14.98 

0.23 


Table 4.3: Results for critical nucleii of Water ( S LW = 15dynes/cm, S p — -50 
dynes/cm ) 


El! 

Po 

r c 

(A) 

AG c l a ss 

(ergs /cm) 

K 1 

(A) 

h 0 1 

(A) 

AG critical 
(ergs / cm) 

@max 

(degrees) 

A U~ t 

&Gclasa 

1 .03 

182.47 

0.846 x 10~ 4 

2.190 

95.840 

0.796 x 10" 4 

52.05 

0.94 

1.10 

56.59 

0.261 x 10" 4 

2.308 

31.012 

0.201 x 10~ 4 

39.08 

0.77 

1.15 

38.59 

0.177 x 10" 4 

2.406 

21.092 

0.109 x lO" 4 

31.72 

0.61 

1.20 

29.58 

0.135 x 10- 4 

2.524 

15.592 

0.610 x 10~ 5 

25.11 

0.45 

1 .25 

24.17 

0.108 x lO " 4 

2.669 

11.870 

0.324 x lO" 5 

18.90 

0.30 

1.30 

20.56 

0.898 x 10~ 5 

2.865 

9.015 

0.146 x 10~ 5 

12.81 

0.16 

1 .34 

18.43 

0.777 x 10“ 5 

3.097 

7.014 

0.557 x 10~ 6 

7.72 

0.07 

1 .35 

17.97 

0.748 x 10" 5 

3.177 

6.515 

0.389 x lO" 6 

6.36 

0.052 

1.38 

.. 

16.75 

0.616 x. 10” 5 

3.604 

1 4.742 

0.305 x 10~ 7 

1.50 

0.0047 
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